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Abstract 

We define a new variable flavour number scheme for use in deep inelastic scattering, motivated by 
the need to consistently implement high energy resummations alongside a fixed order QCD expansion. 
We define the DIS(x) scheme at fixed order, and show how to obtain the small x coefficient functions 
and heavy flavour matrix elements to leading order in the high energy resummation. We then implement 
these results in a global flt at LO which includes leading resummations with running coupling corrections. 
Finally, we address the impact of the resummed results on predictions for the longitudinal structure 
function. We flnd that they stabilise the behaviour of Fl at small x. Overall, we flnd that resummations 
significantly improve the fit to scattering data in the low x regime, although higher orders in the fixed 
order expansion are needed to describe current structure function and related data over the complete x 
range. 



1 Introduction 

The high partonic centre of mass energies available at current and forthcoming particle collider experiments 
impart a greater phenomenological significance to the production of heavy quark flavours. In deep inelastic 
scattering (DIS), for example, heavy flavours can be produced by boson gluon fusion pP once the partonic 
centre of mass energy reaches 4M^, where M is the mass of the heavy quark. For momentum transfers 
Q satisfying 3> M^, one may take the heavy flavours into account by defining heavy parton distributions 
obeying the massless evolution equations When is comparable with M^, however, a consistent treat- 
ment of threshold effects must be implemented. Such descriptions have been around for a long time as so 
called general mass variable flavour number schemes [3 01 |H1 [71 |H| , and have been refined alongside the 
onset of data on the proton structure functions at higher [S] ^1 ^2 ^1 113L I14L [T5] . 



Another problem encountered at high energy is that the coefficient functions relating the proton structure 
functions to the parton distribution functions depend on the logarithm of the Bjorken x variable. So also do 
the splitting functions which govern the evolution of the partons. At sufficiently low x it may be necessary 
to resum log 1 /x terms, thus supplementing the traditional perturbation series ordered in fixed powers of as- 
There is already evidence that such a resummation is necessary in order to improve QCD fits to scattering 
data at fixed order. For example, NNLO fits seem to be improved by the addition of higher order terms 
involving phenomenological small x logarithms |16j . whose coefficients are determined by the data. The 
longitudinal structure function obtained from the reduced cross-section measured at HERA appears to be 
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Figure 1: 0{as) contribution to the boson gluon fusion process whereby a heavy quark pair is produced 
in the DIS final state. 



inconsistent with the theoretical prediction using NLO QCD at small x ^7], indicating the importance of 
higher order contributions. Results for the MS scheme three-loop coefficient functions for F2 JHI and Fl 
|19j demonstrate an apparent perturbative instability in both quantities at low x. In order to resum the 
high energy contributions one uses the BFKL equation |20j . an integral equation for the unintegrated gluon 
4-point function whose kernel is presently known to NLL order |21l I22| . 

The aim of this paper is to present a scheme for the description of heavy quark flavours in DIS, applicable 
to both the small x and fixed order expansions. We implement our scheme in a global fit to scattering data 
at LO which, in the description of F2, incorporates small x resummations at leading logarithmic (LL) order 
with running coupling corrections. We also present resummed predictions for the longitudinal structure 
function Fl, which at low a; is a sensitive probe of the gluon distribution. At high inelasticity y, one cannot 
measure F2 directly but instead obtains a combination of F2 and Fl (the reduced cross-section). It is there- 
fore necessary to have good theoretical predictions for the longitudinal structure function. Furthermore, it 
is widely hoped that measurements of Fl will discriminate between different low x models |23l II 7| . 

The structure of the paper is as follows. In section 2 we briefly summarise variable flavour number schemes 
to introduce our notation, before outlining our proposed scheme (the DIS(x) scheme) in the fixed order 
expansion and at small x to LL order. In section 3 we present our implementation of the small x resummation 
adopting the method of |24| . which includes running coupling corrections. In section 4 we present the results 
of a global fit at LO which includes the resummed heavy flavour contributions. We discuss how this differs 
from a NLO QCD fit. Finally, in section 5, we discuss our resummed predictions for the longitudinal structure 
function. 

2 The DIS(x) Scheme 

2.1 Variable Flavour Number Schemes 

Heavy quark pairs in DIS can be produced by boson gluon fusion T - the leading order diagram at the 
parton level is shown in figure H"^- Higher order contributions can involve the gluon interacting further with 
quarks from the proton, as well as real and virtual emissions. Taking into account all such diagrams, the 
heavy flavour contribution to the proton structure functions {Fa} can be written as: 

rif 

P-'" = ^ E + C£f,) ® S + Ci^^ ® g + nfCZI, ® Qns] + ® E + Cf,^, ® g] , (1) 

f i=l 

^Here we are concerned with heavy quarks which are generated perturbatively, as opposed to intrinsic heavy flavour distri- 
butions, which would be present in the proton at all values of (Sjj- 
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Table 1: The leading order of each of the coefficients in the heavy flavour structure functions. 



where n/ is the number of active light flavours (n/ = 3 in phenomenological applications). Following p], 
we have explicitly decomposed all coefficients into singlet and non-singlet parts, and also separated the 
coefficients into parts arising from the photon coupling to a light quark pair {Ca.qi{x, Q"^ /M'^)) from those 
where it couples to the heavy quark pair ^ {Ca.Hi{x,Q'^ /M'^)), i € {q,g}- The coefficient Ca,Hq has only a 
pure singlet contribution, as only an intermediate gluon can couple the heavy quark pair (which interacts 
with the photon) to a light quark pair. This is a flavour singlet exchange. The leading order of each coefficient 
is given in table ^ One could choose to use this description of the structure functions within a particular 
factorisation scheme at all values of Q^. This is the fixed flavour number scheme. Collecting together the 
coefficients of each parton density, one may tidy up equation 



P _ ^FF 



s + c 



^^(FF)T, 
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defining the fixed flavour coefficients f and 
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c. 



c. 
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However, a problem arises at high in that the fixed flavour coefficients are logarithmically divergent in 
Q'^/M^. For example, one has in the MS scheme: 



^FF TP 



(4 - 8a; + Sx^) ( log + log(l - x) - log a; ) - 4 + 32a; - 32x'^ 



+ Oial), (5) 



with Tn = 1/2 a colour factor. In principle both the partons and the coefficient functions (beyond leading 
order in as) are dependent on the coUinear factorisation scale We adopt the natural scale choice 

/^F ~ from now on. The problem of divergent FF coefficients can be circumvented by defining a heavy 
quark distribution at a matching scale — appropriate to the threshold for heavy quark production. 
Singularities are then absorbed into the parton distribution, and its evolution is governed by a splitting 
function as in the DGLAP equations for light quarks. Below the matching scale, there is no heavy quark 
distribution. Different choices can be made for the matching scale - we use = Af^. This ensures that up 
to NLO in the fixed order expansion, the heavy partons evolve from a zero initial value. This choice also 
reduces technical complication at higher orders, although no choice of matching scale will ensure continuity 
of all partons. One now has a variable flavour rmiBher scheme, as (the number of active quarks) increases 
as crosses the matching scale from below. Equation Q applies for < A-P , whereas for > one 

^There are technical complications arising at NNLO due to heavy quark pairs produced away from the photon vertex, which 
contribute to both the light and heavy structure functions These do not concern us in this paper. 
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has: 

= C(^^)T®/("^+i\ (6) 

where qn (qh) is the quark (antiquark) distribution for the heavy flavour. Such a scheme was first suggested 
in 0] , and a proof that the variable flavour number scheme is well-defined to all orders in perturbation theory 
was given in j27| . 

The general requirements a variable flavour scheme must satisfy were derived in [S]. The fact that the heavy 
flavour structure functions cannot depend upon whether one uses a fixed or variable flavour number scheme 
imposes a relationship between the variable and fixed flavour partons. For the heavy flavour and singlet 
densities: 

= i 1 = 1 A,, a;, k x:; )=a/("/), (7) 




where the heavy flavour matrix elements {^y} are perturbatively calculable in a given factorisation scheme, 
and qH+ = Qh + qn- They satisfy: 

A, =S,,S{l-x)+Oias) (8) 

n.b. at leading order, the (jif + l)-flavour parton distributions are the same as their n/-flavour counterparts. 
Similarly, non-singlet quark combinations are related by: 

Juf + l) ^ ^JVS ("/) .g^ 

As Q^/M^ oo, the heavy flavour mass becomes negligible in relation to the scattering scale, and so 
we must be able to neglect terms ^ 0{M'^ /Q'^). Thus, the variable flavour coefhcients must tend to the 
appropriate massless coefficients in the asymptotic limit of high Q^. Like their fixed flavour counterparts, 
they are only defined up to a factorisation scheme choice. Equivalently, one is free to transform the singlet 
parton densities according to: 

where the matrix Z in each case has diagonal elements d{l — x) + 0(as), and off-diagonal elements 0{as)- 
There is a further constraint from the momentum sum rule: 

/■I 
Jo 

where T,(x) is the quark singlet distribution: 

Eix)^J2^q,ix) + q,{x)]. (12) 

i 

The sum is over all active flavours. A similar transformation to l)f 0|1 also holds in the non-singlet sector, but 
with a 1 X 1 transformation matrix. 

There is a further ambiguity in the VF coefficients arising from the fact that there is one more VF coefficient 
than there are FF coefficients. Combining equations (001 one finds: 

Fa,H = iCl% ® Auq + Cl^-''^ ® A^^ + Ca,g ® Agg) ® 

+ {Clfj ® Ang + Cl^^'''' A,g + Cl^ ® Agg) ® gM, (f3) 
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Equating the gluon coefficients in equations ((Tl 



FF 



^a,H ® + C^,q'^^ ® Agg + C'^j ® Agg — C^ g . 

For example, using the fact that Agg = 5{1 — x) at leading order, one obtains at 0{as)'- 



C. 



FFil) 
2.9 ' 



where in the MS scheme: 



4(1) _ T 



(4-8x + 8a;2)log 



91 

AP 



(14) 



(15) 



(16) 



is the coefficient of as / (47r) . Comparing with equation |SJl , one sees that the heavy flavour matrix element 
contains the same divergence as the fixed flavour coefficient, so that the variable flavour coefficient is finite 
as — > 00. 

Equation p4(l shows that although each VF coefficient must tend to an appropriate massless limit, one can 
choose to shift terms vanishing at high from into C^J. This can be examined in more detail in 

Mellin space, defining the moments of parton distributions ^ by: 



q{N) = MN[q{x)] = / dxx'^qix) 



(17) 



Heavy flavour coefficients contain an implicit Heaviside function Q{W^ — AP), where is the proton- 
photon centre of mass energy (n.b. must exceed the threshold for pair production). Correspondingly, 
there is a maximum value of x: 

4Af2^ 



1 + 

We thus define moments of the heavy coefficients via the scaled variable x' ~ x/xmax to obtain: 

MN[f{x',Q^/M^)] = C dx'x^'fix'^Q^/NP). 
Jo 

One then has the convolution theorem: 

MN[f{x',QyAP)<g)g{x,Q^)]=MN\ f f 



(18) 



(19) 



x' 



giz,Q )dz 



z ' M2 

N[f{x',Q^lAP)]MN[9{x)l 



(20) 



where N is conjugate to x' on the left-hand side, but on the right-hand side is conjugate to x' in the first term 
of the product, and conjugate to x in the second term. Taking x' moments of equation (|14|l and rearranging, 
one obtains the A^-space gluon coefficient: 



r^VF _ ^ 

^a.cj - J- 



/=.FF /^VF,PS A r<v J- A 



WF 



(21) 



Substitution of the perturbative expansions for the coefficients and heavy flavour matrix elements leads to a 
power series in as, which can be inverse transformed term by term to give the resulting a:'-space expression. 
Choosing the quark coefficients then fixes the gluon coefficient according to equation (|21|l . 

To summarise, there are two possible types of choice in specifying a variable flavour number scheme for 

^Another definition of the Mellin transform is commonly used, corresponding to Mjv-i in our notation. 
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Figure 2: Diagrams contributing to the LL impact factor. Two further diagrams are obtained by reversing 
the direction of the quark loop. 



1. Method of colhnear factorisation. Moving between different choices involves an explicit transformation 
of the partons according to equation H10() . Such an ambiguity is always present in the coefficient 
functions, for both the light and heavy flavours. 

2. Position of M"^ /Q"^ terms. Choosing a different definition of e.g. the gluon coefficient does not involve 
a transformation of the parton densities. 

Consider a change of definition of type (2) above, within a fixed colhnear factorisation scheme, from a scheme 
(A) to another scheme (B). Up to ©(al), one may ignore the terms involving the pure singlet coefficient 
and solve explicitly from equation (|21|l for the A^-space gluon coefficient in the new scheme: 

c^^r' - [^IT^ - c^a!'/']^ + ciT'^ (22) 

^gg 

where the bracketed term is ©(M^/Q^), and to leading order one may set the denominator to 1 on the 
right-hand side. 



2.2 A VF Number Scheme for Small x Physics 

We now turn to the problem of how to formulate a VF scheme that can describe heavy quark production at 
small X, but also matches on to the fixed order expansion. One must first satisfy requirements (1) and (2) 
above, namely choose a suitable factorisation scheme and positioning of /Q^ terms. At small x the heavy 
flavour contributions to the proton structure functions are given by the high energy factorisation formula 

EHll2n|: 

Fa,H^asJ^ ^K{eiQ\Q^lM^)f{N,e,Ql)~gB{N,Ql), (23) 

where /(A^, fc^, Qq) the unintegrated gluon density, and gsiN, Qq) a bare gluon distribution which absorbs 
colhnear singularities upon solution of the BFKL equation. The {ha} are impact factors coupling the virtual 
photon to the gluon via a heavy quark pair (see figure whose LL forms may be found in |30[ I29| . The 
corresponding results for the light structure functions are in j31| . The NLL results are not yet known, but 
work is in progress [32 ESI ElESl EOl E3 • 

It is tempting to immediately identify the impact factors with the heavy flavour coefficient functions C2.Hg 
and CL,Hg, both below and above the matching point = M^. However, whilst the longitudinal impact 
factor is finite as Q^/A'P — > oo, the factor h2 diverges as a consequence of the fact that F2 is proportional 
to the quark singlet distribution which is intrinsically non-perturbative {Fl, on the other hand, begins at 
0{as))- Instead one may consider the quantity: 

-as %h2{kyQ^QyM^)me,Ql)9B{N,Ql), (24) 



dlogQ"^ Jq k 
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which serves to define the factor /i2. The resulting coefficient is then not simply C2.Hgi but in general is a 
mixture of the coefficient function and anomalous dimensions 7_ffgj 7.33 entering the DGLAP equation and 
governing the evolution of the charm and singlet distributions. The interpretation of this quantity depends 
on the factorisation scheme. In the massless case (recovered here by taking Q^/Af^ — > 00), it is simplest to 
adopt the DIS scheme in which the structure function F2 has the naive parton model form: 

"/ 

F2(x, Q2) ^ Y}^}[q,{x, Q2) + Q^)] (25) 

i=l 

to all orders, where ei is the electromagnetic charge of quark species qi. The impact factor h2 can then be 
interpreted as the quark-gluon anomalous dimension 7^^^^. If one were to adopt this in the massive case, 
however, the massive impact factor /12 would be interpreted as an M^-dependent anomalous dimension cou- 
pling the heavy flavour distribution to the gluon. This is undesirable for a number of reasons. Firstly, there 
is difficulty conserving momentum either side of the matching point if splitting functions of massless quarks 
differ from those of the heavy species. Furthermore, mass-dependent splitting functions will lend consider- 
able analytic complication to the fixed order expansion. The massless MS scheme anomalous dimensions 
are well known (now up to 0(0^5) ISHEHDj and can be straightforwardly transformed to other schemes. 
More importantly, the use of massive splitting functions for the heavy partons makes them different to the 
massless partons. Consider the structure functions at high . The fact that the massive partons have 
evolved differently to the massless ones from their matching scales up to high means that the coefficient 
functions are also different, and do not reduce to the conventional massless coefficients Such a division 
between the parton species seems unphysical. 

Instead it makes sense to factorise the partons in the DIS scheme such that both the heavy and the light 
ffavours evolve according to the massless DIS splitting functions. One can ensure this by transforming from 
the MS scheme partons. The explicit iV-space transformation for < Af ^ is entirely mass-independent and 
takes the form (see e.g. |42) ) : 

where ffavour factors have been explicitly displayed in the massless MS scheme coefficients. This leads to a 
similar transformation matrix for the singlet and heavy quark distributions above the matching point^: 




2inf + l)C^^ C^^ 



( 

5^ |, (27) 



with qH+ = Qh + Qh ■ For the non-singlet sector: 



Qris — qns ■ [^o) 

The DIS scheme partons now have no mass dependence, other than that introduced by the transition points 
— AP for each heavy flavour. Instead, explicit dependence on M'^/Q'^ resides in the FF and VF coeffi- 
cient functions, rather than the anomalous dimensions governing the parton evolution. The transformations 
above explicitly define a coUinear factorisation scheme, thus fulfilling condition (1) above. 

One must now specify the position of 0{K'P /Q^) terms. It has been suggested previously in the MS scheme 
fixed order expansion .43. that one should choose the leading order contribution to the heavy quark 

^Strictly speaking there is a further pure singlet contribution involved in the transformation of the single quark species qH+- 
This can be neglected up to ©(og). 
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coefRcient to be C. 



MS(0)yF 
2.H 



5{1 — x'). This is known as the ACOT(x) scheme. Note that this automatically 

c. One can extend this to all orders i.e. 



generates the correct asymptotic limit at high ^ where x' 
^2^h'^^ — C^^'qW) to give a simple, physically motivated definition of the quark coefficient at all orders. A 
similar definition can also be chosen for the heavy quark coefficient entering F^^h |44j . 



^DIS.VF 



We propose using the same specification for the DIS scheme. Here, this amounts to choosing ^ — 
6{1 — x') to all orders in as- It will be seen that this simplifies considerably the solution of the gluon 
coefficient in the small x limit. This is clearly a physical way of defining the heavy quark coefficient within 
the DIS scheme, which is itself a more physical scheme to work in due to its connection with the naive parton 



DIS V F 

model. The natural choice for F^^h is also to choose ^' 



CY^^{x'), where the right-hand side denotes 



the functional form of the massless coefficient. However, a problem occurs with smoothness of Fl^h across 
the matching point. At 0{as) in the longitudinal case, equation H14|) gives: 



^VF{1) 



= C 



FF{1) 
L,g 



and one has pQ: 



with: 



C 



FF(1) 



4x(l 



2v-{l- u2)ln 



1 - i; 



e{W'^-AM'^), 



UPx 
QHl-x 



(29) 



(30) 



(31) 



The parameter v is the velocity of the heavy quark or antiquark in the photon-proton centre of mass frame, 
and hence satisfies < v < 1. Taylor expanding H3U|I gives: 



C 



FF(1) 



8x{l~ x)v^ + Oiv^). 



The fixed and variable flavour descriptions of Fl h dictate at 0{as)- 



F, 



FF 
L,H 



(^)<7)«,(-)+o(4); 



pVF 



\ 4.11 J 



C 



a 

VF(1) 

i:3 



I + '^L.H 



' qH+ 



0{a, 



s)- 



(32) 

(33) 
(34) 



The FF contribution and the first term of equation (|34|l are suppressed by , whereas the second term in 
equation l|34(l is not if one chooses C]f^ — CY^q{x'). Even though the heavy quark distribution will be small 
at the matching scale, one will still have a kink in the structure function at = ^vf^. Smoothness of the 
structure function can be recovered by instead demanding C]^^ = f{Q'^/AP)C^^^{x'), where f{Q'^/M'^) is 
a function such that /(I) — and /(oo) = 1 with a smooth interpolation between these values. A suitable 
choice is |44| : 

/Q^\ 5/ 1 1^ 



/ 



1 + 4Af2/Q2 



(35) 



which fulfils the above conditions and is suppressed by the heavy mass factor 4Af2 at intermediate values. 

is not affected by this choice, but is the same as that obtained by choosing / = 1. This 



Note that C 



VF{1) 

L..g 



(B)(1) 

i:3 



when transforming between schemes dif- 



foUows from equation (j^ . which at 0{as) gives C)^ g = C 

fering only in the placement of 0{M'^ /Q"^) terms. The NLO coefficient C"]^^'^' and higher orders, however 

will be affected as they depend upon the choice of C]^^^^^ ■ Note the prefactor of equation H35|) is chosen to 
have no dependence on x, thus introduces no extra complication in convolutions involving the heavy quark 
coefficient, matrix elements and the partons. 



There is a further problem in -Fb,//, regarding the order of the perturbation expansion cither side of the 
matching scale. For < M^, one has F2,h 0{as), whereas i^2,ff C>{ag) for > A/2. One must 
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choose how to deal with the ordering in such as way as to guarantee continuity of F2,h in the fixed order 
expansion. FoUowing we add the frozen contribution: 



P2,H°^ = ^2%(^, = M') ® 9{^, = M^) (36) 
for > (see |33 for a discussion). 

One has now satisfied choice (2) above for both F2^h and Fl^h, and thus completely defined a variable 
flavour number scheme. By analogy with ACOT's proposal, we refer to this as the DIS(x) scheme. In the 
next section, we will present the derivation of the small x variable flavour (and fixed flavour) gluon coefficient 
functions. To use them alongside the fixed order expansion, however, the known MS scheme heavy flavour 
coefficient functions and matrix elements must be transformed into the DIS scheme. 

Transformation of the FF coefficients proceeds from equations 1261 1^5)) . together with scheme independence 
of the structure functions: 

FaM = C'^^'^IS) T ® /("/)DIS ^ fj^FFMS) T ^ f{nj)MS^ (37) 

One obtains: 

( C^^^^^ ) = ( Cl^MS c^FME ) ^ ( z^^"^^ ^ ^3g^ 



2,ns 



where Z^"^-* is the matrix of equation l|26|l . Here inversion of the matrix under a convolution sign can be 
understood by transforming to iV-space, substituting in the perturbative expansions of all quantities, then 
inverse transforming term by term back to x'-space. A similar procedure can be carried out for the variable 
flavour coefficient functions, and one obtains: 

( C^;^'™^ (jVF.ms (jVF,ms -j ^ (■ C^^^ (jVf.ms (jvf.ms )[z("-f+^)]"\ (39) 
with Z^"^"*"^-* the matrix of equation H27|) . Explicit results up to 0{as) for the FF coefficients are: 

^F^^,DIS(1)^^™(1). (40) 

^F_F,DIS(1) ^ ^FF,MS(1) ^^^^ 

For the VF coefficients: 

C™^'^ = 0; ' (43) 

^yF,DIS(l) _ ^VF,MS(l) r,^MS(l) 

^2,g — ^2,g ~ '^^2,g i [^^) 

^VF,DIS(1) ^ ^VF,MS(1). /^gN 

^VF.msii) ^ ^yF,Dis(i) (4g) 

Note that to this order one has 6*2^^'°^^ = G^!h^x') as required by the choice of where to put 0{KP/Q'^) 
terms, thus there is no need to perform a further transformation according to equation (|22|) . 

One must also transform the heavy flavour matrix elements {^ij} from the MS scheme, where they are 
known up to 0{a'g) to the DIS scheme. This can be done using scheme independence of F2, as well as 
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equivalence of the FF and VF descriptions: 

^ ^(VF,DIS) T ^ y^DIS ^ ^(n/)DIS 



and therefore: 

Exphcit results up to 0{as) are: 



^Dis ^ <^ ^ (47) 

(48) 

A°fW=A^(i)-2Cy«; (49) 
<r=^S^^^+2Cr^'. (50) 



Expanding equation (|14|l to 0{as) in the MS scheme, one finds: 



^2,g — '-'2,9 ~ ^ifg ■ \^ ) 

Applying the results of equations BUI 1501) one then finds: 

^FF,DIS(l) _ ^yF,DIS(l) 4DIS(1) /,r,N 

which is exactly what one expects in the DIS scheme, given that in equation H14|l no factorisation scheme is 
specified. This serves as a consistency check of the above transformations. 



2.3 The DIS(x) Scheme at Small x 

To see how the above scheme is implemented in the small x limit, it is convenient to Mellin transform 
equation according to: 

F{-f,N)= dk^{k^)-^-^'F{k^,N). (53) 



We now drop tildes in N and 7 space, instead denoting explicitly the arguments of each function. For F2^h'- 

F2,Hh, N, Q^/M^) = /i2(7, N, Q^/M^)f{N, 7, QDgsiQl N) 

^h2{i,N,Q^/M^)g{N,^). (54) 

For < M^, one identifies: 

Cl^{l, N, QVM^) = h^ii, N, Q^/M^). (55) 

For > M^, one must consider the VF expression for F2 in the DIS(x) scheme. As has already been noted, 
the impact factor /12 diverges as Q^/M^ 00 due to lack of coUinear safety, thus one cannot directly take 
the variable flavour coefficient from /i2. However, one can make progress from equation (|13|l expressing the 
relationship between the VF and FF schemes. Ignoring the pure singlet term in equation (|14|l . and using 
Agg = (5(1 — x) and Cj'^ = (5(1 — a;') at LL order (where the second equality follows from the definition of 
the DIS(x) scheme): 

Angii, N, Q^M') = CiJ{^, N, Q^M^) - CYJ{^, N, Q^M^) (56) 
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Given that the cquahty holds for any value of Q^, one may consider the limit 



One must have: 



0, 



(57) 



i.e. the variable flavour gluon coefficient tends to the appropriate massless limit at high Q^, which is zero. 
Then the second term in equation (|56|l vanishes and leaves: 



AHg(7, Q'/M') = Cl'g (7, TV, Q'IM')\ 

where the notation on the right-hand side implies that the parts which vanish as Q^/M^ 
Thus, the DIS(x) VF gluon coefficient at small x is given in double Mellin space by: 



(58) 

oo are removed. 

(59) 



It is not immediately obvious that all of the requirements of the DIS(x) scheme have been satisfied. It 
must be checked that the partons evolve according to the massless DIS scheme splitting functions, and 
that the heavy quark coefficient is trivially defined to all orders. One way to show that these conditions are 
implemented successfully is to derive the above result in a different way, which highlights their role explicitly. 

Considering the derivative of F2.H with respect to logQ^, one has: 



5F2,h(7,^,QVM') _ dCl^{j,N,Q'/M') 



91ogQ2 \^ 91ogQ2 , v,,^-,^ ^ -2,«V,,^-,^ 

+2inf + l)<f ''^(7, N, gVM')7,s(7, N, Q^/M^)) 3(7, iV) + . . . , 
^h2[l.N,Q^/M^)g{-1,N), (60) 
where the ehipsis denotes terms in partons other than the gluon. The final term in the bracket can be 

PS V F 

neglected to LL order, as at lowest order in the small x expansion "fqg and ^' both contain NLL terms. 
In iV-space, these correspond to contributions of the form ^ as{cy-s / N)'"" . The product of two such terms is: 



as 



\N J 



2 



which gives a NNLL contribution. The impact factor h2 is related to h2 by: 

dh2{j,N, qyAP) 



h2{i, N, Q'lM') = 7/12(7, N, Q'/RP) + 



9iogg2 



(61) 



(62) 



One now implements the DIS(x) scheme by identifying '^qg{'^, N, Q'^/AP) with the limit of /i2(7, N, Q"^ /Iv'P) 
as /AP — ^ 00 (this corresponds to the massless DIS scheme anomalous dimension 7^^^(7, N)), and setting 
C2Hi'y,N,Q'^/AP) = 1. Rearranging equation (|60|l then gives: 



dCYJ{j,N,QyAP) 
d log Q"^ 



+ iCI^[i.n,q^/aP) = 



9/12(7, N, Q^AP) dh2{i, N, Q^/AP) 



+ 7 



9iogg2 aiogQ2 

/i2(7, Q'IAP) - /i2(7, A^, QVm')I 



(63) 



to be solved subject to the boundary condition of equation (|57|l . The solution is then given by equation 
(I59|l . as found previously by considering the equivalence of the FF and VF descriptions for ¥2- The former 
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method exphcitly demonstrates continuity of the structure function at small x across the matching scale 

= M^. The second method involving the derivative of F2 allows one to see how the definition of 
relies upon the DIS(x) scheme choice of splitting functions and heavy quark coefhcient. 

For the heavy flavour contribution to the longitudinal structure function, one has in double Mellin space for 

FL.H^hL(-i,N,Q''/M'')g{-f,N). (64) 

For > M^, implementation of the DIS(x) scheme is somewhat simpler than in the case of i^2,_ff- The 
longitudinal impact factor is not divergent as Q'^/M^ 00, and so one can apply the factorisation formula 
(jMjl at aU scales. Thus the coefficient C^J(7, A^, QV^^) at small X is obtained from the impact factor 
/il(7,^,QVm2) asforCf^f. 

There is potentially a discontinuity in Fl^h across the matching point in the small x limit, arising from the 
non-zero initial value of the heavy parton distribution. Analogously to equations H33I34|I one has: 

Ff_|;(x', asiQ'), QVm') = C^ix', Q^/M^) ® g{x, Q'); (65) 

FlM^', o^sm, QVM^) = CYUx', QVaP) ® qH+ + Cll(x\ qyM^) ® g{x, Q^) 

= CY'nix', QVM^) <g> Angix', Q^M^) <g> g{x, Q^) + C^^(x', Q^M^) (g> g{x, Q^). 

(66) 

At = A'P , the FF contribution and the second term in equation (|66|l are equal, as the FF and VF gluon 
coefficients are both obtained from the same impact factor. The potential discontinuity arises from the first 
term of equation H66|l . although this is formally 0(a|(as/A^)") in the high energy expansion given that Ahq 
and Cl,h both contain terms in as{as 1^) as their leading high energy divergence ^. However, the prefactor 
/(Q^/Af^) of equation H35I) suppresses the first term of equation H66|l . which is then zero at = itself, 
thus facilitating smoothness of the structure function. Nevertheless, continuity of F^jj is only formally 
restored by taking higher orders in the resummed coefficient. 

Care must be taken in implementing the resummed VF coefficient, due to the frozen term at — M"^ 
used in the fixed order expansion (equation H36|l 'l. The VF coefficient implicitly contains the fixed flavour 
coefficient according to equation l|59|l . Thus, one must subtract the 0{as) contribution to the resummed 
^2g from the resummed VF coefficient to avoid counting this term twice at = M^, as it is already 
included in the fixed order expansion. Explicitly, the FF description at LO and LL order gives: 

fIh = (x',^) «5("^^(x,Q^)+C,^f- ®9'-H-.Q'\ (67) 

where: 

is the total resummed FF coefficient with its leading term subtracted (n.b. this remainder is 0(a|)), and 
C'2g'^^'^^^ denotes the 0{as) contribution to the resummed coefficient. One subtracts the ©(as) term from 
the resummation, as this is already included in C2 ^ . Naively in the VF description, one would write: 

^g("/+i)+C2';f"®<?("^+i\ (69) 

similar problem occurs in the fixed order expansion, where there is a discontinuity at 0{ag), due to the non-zero value 
of the NNLO matrix element A^)^ at = (n.b. at 0(a^g)). See for a discussion. 



pVF 
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where the frozen term is included in the fixed order expansion at 0{as)- At LO and LL order, the («/ + 1)- 
flavour partons are the same as their n^-flavour counterparts, and using equation (|57|l in x-space, one 
obtains: 



VF 



2,H 



A 



Hg 



as 
47r 



FF{\)i , 
2,3 



FF,rs{\) ( I 



2,3 



"9 I ^ ' M2 



FF,rs 
2,3 



M2 



(70) 



where the first hne originates from the heavy quark. Evaluating this at = M"^, one obtains: 



_ fas 



2.H 



L 47r 



<f 1) + C.^f'^^^^^lx', 1)] + Clf^\x', 1)} ® M\ 



(71) 



Comparing this with equation , one sees there is a double-counting of the high energy limit of the LO 
FF coefficient, present in both the frozen term and the VF resummation. One must therefore subtract 
C2g''^'^^^\x' , 1) (where x' is evaluated at = M^) from the resummed VF coefficient. Equation then 
becomes: 



as 
47r 



VF.i 
2,3 



,9 

,("/+!) 



M2 



'5 



("/ + !) 



(72) 



2.4 The Heavy Flavour Matrix Element Aug 

A further check of the above results is obtained by considering the LL limit of Aug , as given by equation 
(|58|l . Using the form of the impact factor /i2(7, N) [30], one finds: 



Aug (7, Q^/M^) = h2{7, N, Q^/M^)\ ^2 
3/1(7) 41-^ 



272-7-1/02 



167r2 (7-57)(l + 27) 



V 4Af2 



V^(2 + 37 - 37^) r(7) 



r(i/2 + 7)j ' 

(73) 



with: 



4^ 7-57 r3(i-7)r(i + 7) 

^^^^ = T"^3-27 r(2-27) ■ 



Expanding equation (|73|l in 7 yields: 



(74) 



AHg{i,N,Q'/M')='^ 



4 g2 2 

3^°Sm^+3 



, 2 2 20 Q2 

7 (-3 log ]J7I + ylog]^^-^ 



10 27r2 
9~ 



0(7^) 



At LL order with fixed coupling as, 7 is identified with the BFKL anomalous dimension |20j : 



(75) 



as , 

7BFKL — + 



2C(3) (^^ 



TV / 



(76) 



and so the high energy N-space heavy flavour matrix element is given by: 



AHgiN,Q'/A\Q'/M') 



as 



4 g2 2 

3^°Sm^ + 3 



/asy 
V47r/ 



2 Q2 ^ 80 g2 4Q 8 ^ 



TV 



+o(4). 

(77) 
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This can be compared with the N ^ hmit of the fixed order matrix element in the DIS scheme. For the 
0(0^) contribution, one needs the result (derived in x-space from equation 1471) '): 



"2.9 



A 



MS(1) 



(78) 



The corresponding result for ^ is given in equation (|50|) . At 0(a|), A||g can be found in jSj. Applying 
the above transformations yields: 



2 

M2 



= 2 
_ 1 

X 



A[x^ + (1 - a;)2] loj 



01 

Af2 



-81og^ 



Af2 



80^ Q2 



(1 - 2.x + 2a;2) log 
40 Stt^ 



1 - X 



- 1 + 8a;(l - x) 



+ 0(x°). 



(79) 
(80) 



Taking the Mellin transforms M[Ar[A^^ ] then the limit ^ reproduces the terms in equation H77|) . 
Thus the LL terms in the heavy flavour matrix element are correctly predicted by the small x resumniation. 
It is instructive to examine the behaviour of A^g when — j\,f2^ ^j^-g governs the initial value of the 
charm distribution at the matching point. It may be seen that the term in in equation (|75l) is negative at 
(32 = potentially dr iving the entire matrix element negative as A^ — *■ resulting in a charm distribution 
which does not turn on with a positive value. This is not formally a problem given the interpretation of 
parton distributions as probability densities is only true in the naive parton model. Such a situation already 
occurs (and can be dealt with jHlEI) in NNLO perturbation theory in the MS scheme, where the 0{asY 
contribution to Aug 5 ensures that the charm turns on with a negative value. However, one would hope in 
a more physically motivated scheme such as the DIS(x) scheme that the partons are positive, due to their 
being related at zeroth order to the heavy flavour structure function F2.h- 



In fact, one may resolve this by considering higher orders in the fixed order expansion at small x. Evaluating 
the numerical coefficients of the leading logarithmic terms yields: 



A 



Hg\Q2=M2 = a 



.66667 - 12.986-^ + 267.96 ^ - 2794.8 ^ + 31076 
+1946155 (^yVc'(a' 



(81) 



where a = as/(47r). One sees that not every term in the fixed order expansion is negative. At ©(q?!), 
the complete matrix element becomes positive again in the fixed order expansion at small x. The LL or- 
der resummed matrix element expanded to 0{a^g) (a very good approximation to the all orders result for 
x > 10"^) at (52 — for M — 1.5GeV is shown in figure |2| Although slightly negative at higher x, the 
small X limit is positive and sets in at even intermediate values of x. One must also bear in mind the 6{l — x) 
contribution, absent in the plot, which is also positive. Indeed, the matrix element is positive for all values 
of A^ in Mellin space. The positivity as a; ^ is in contrast to the conclusion reached at O^Og) in the 
fixed order expansion. This highlights the importance of the high energy resummation, and it is reassur- 
ing that such a physical scheme choice results in a positive heavy quark contribution to the structure function. 



3 Resummed Heavy Flavour Coefficients 
3.1 The BFKL Equation with Running Couphng 

So far we have presented the DIS(x) scheme at fixed order, and in double Mellin space at small x. Explicit 
results for the x'-space heavy flavour coefficients can only be obtained after solving the BFKL equation for 
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the gluon density in equation (|23|l . In this section we adopt the method of solution of which is at LL 
order but includes the running of the coupling as- Although formally a NLL effect, including the running 
coupling changes the nature of the solution of the BFKL equation and helps to stabilise the small x expansion. 

A comment is in order regarding the use of resummed coefficient functions and quark-gluon splitting func- 
tions in a LL fit. The impact factors calculated at LL order actually contribute to the coefficient and splitting 
functions at NLL level (n.b. terms of order as{as/Ny^ in Mellin space), according to the definition of e.g. 
|45| . Only the anomalous dimensions and contain LL terms, which arise from the LL BFKL kernel. 
We adopt the viewpoint that the leading resummations should be present in each physical quantity, and thus 
include the effects of the impact factors in a fit to data. After all, the quark-gluon splitting function and 
coefficient functions are zero in the strict definition of the LL expansion. Not including the NLL terms in 
these is analogous to disregarding the evolution of the partons in a LO fit in fixed order perturbation theory 
due to the fact that only the quark coefficient is non-zero at 0{a^). This is nonsensical if LO QCD is to 
differ from the naive parton model. Likewise, if resummed perturbation theory is to differ from the LO fit, 
in both the parton coefficients and the evolution, one must include the impact factors. 

We now recap the essential features of the method. A full presentation and discussion can be found in j24| . 
When the LO running coupling is used, the BFKL equation in double Mellin space becomes a first order 
differential equation |46j: 

where //(7, Qq) is a non-perturbative input gluon 4-point function, usually taken to be to be a delta function 
in momentum space S{k'^ — Qq). Also, xo is the Lipatov function: 

Xo(7) = 2V'(l)-^(7)-^(l-7)- (83) 

Solving, one obtains: 

/(7, N) = exp [-Xoh)/{PoN)] H ^lltl^ exp [Xo{^)/0oN)]d^, (84) 



7 



with Po = 7r/3o/3, and: 



1/2 



2V'(l)(7-iVln^ ^^^^ 



2j Vr(i-7) 



(85) 



Up to 0{Qq/Q^) corrections, the integral in equation H84() can be taken from to 00 [21], which leads to 
factorisation, up to power corrections, of the integrated gluon distribution as follows. Including the bare 
gluon as before, the integrated gluon structure function is obtained by inverse Mellin transformation of 
equation ifSl) ^: 

G{Q'.N)= 77^- i7 fh,N)gB{QlN). (86) 



l/2-^oo 27ri7 VA2 

The Mellin transform of // = 6{k'^ — Qq) with respect to fc^/A^ is (Qq/A^)^'''"^. Inserting this into equation 
TO yields: 

QiQ^ N) = gEiQ^N)giiQl N)gB{Ql N), (87) 

where : 

1 ^l/2+ioo -. 

g^iQ^N)^— -/'3°exp[7ln(QVA2)-Xo(7)/(/3oA^)]d7 (88) 
27ri^/i/2-zoo 7 

*A factor of 7"^ now arises from considering the integrated gluon structure function. The unintegrated function is related 
to this by a derivative with respect to In corresponding to a factor of 7 in Mellin space. 
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and: 

GiiQl N) - / exp [-7 In {Ql/ + Xo(7)/(^o^)]d7- (89) 
Jo 

Note that a 7-independent factor of — (A^/Qq) In (Q^/A^) has been absorbed into gsiQa, N). The function 
Qi is not calculable due to the presence of branch points on the contour of integration, but evolution of the 
gluon distribution depends only on the perturbative piece Qe- The latter has been modified by a factor 
originating from the choice of renormalisation scale of as ^: 

= exp[l/2(ln(xo(7)) + ^0(7))]- (90) 

To calculate a structure function, one follows the above derivation but includes an impact factor hai^), to 
yield the perturbative contribution (to be combined with the gluon distribution): 

^E,a^^ ' ^/'3n7)exp[7ln(gVA2)-Xo(7)/(^o^)]d7. (91) 

For the longitudinal coefficients and F2,h FF coefficient, one uses the factors Hl and /i2 respectively. For 
the i^2,_ff VF coefficient, one uses the result of equation H59II . Coefficient functions are now obtained via: 

<f'''^^s{Q%N)^^0^. (92) 

The integrals can be calculated numerically, but approximate analytic expressions may be obtained by 
expanding the Lipatov function: 

Xo(7) = - + E2C(2« + l)7'", (93) 

valid for —1 < 7 < 1. One then finds, by deforming the integration contour of equation H91|l to enclose the 
negative real axis, and taking into account the discontinuity across the cuts 24j 



a„7'"+' W7, 



^ 1 ^^^^ 

exp I 7^ - ^ 

Po 

where t = In (Q-^/A^), a„ = 2C,{2n + \)/{2n + 1), and 7^ is Euler's constant. There is a correction of order 
/Q'^ to equation (jHU, given that equation H93|l is not valid for 7 < —1. To evaluate this integral, one first 
changes variables to y = 7i, and then Taylor expands in y the factor: 



^ ^ L n— 1 J n— 



(95) 



Then one has 

Ti? _ = — sin I _ , , _ 



.F^,„ = -sin(_)exp(3-Ui/(/3oAr) 



Y,KaAPo,N)t-- f ^-^/('^"^^-lexp 



(96) 



Each of the integrals in the sum can be evaluated using: 
rO 



dyy-^iM-i exp {y)y- = ^^{-\ / {P^N))T{-1 / {p^N)) , (97) 



'See 1241 for a discussion. 
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where: 

n-l 

A„(-l/(^o^)) = E (-l)"rf™,n(^oiV)-"+'", (98) 

rn—O 

do^n — 1, and dm.n+i — drn.n + ndm-i.n- Setting m — n in this recurrence relation yields dn.n+i = ndn-i,m 
and hence: 

d„_i,„ = (n- 1)!, (99) 

and so the series in equation (|96|l has factorially divergent coefficients Ka,n{^o, at high n. Thus the series 
is asymptotic, and must be truncated at a finite order n = uq. Finally, one has: 

TE,a = - Sin (J^) exp r(-l/(^oiV))tV(^oA^)(_i)-i/(/5oA.)-i 

(100) 

xJ2KaAPo,N) (a5/5o)" An (- 1/ (^oA^)) , 

using = (as/^o)", as = Sas/ir. A similar expression is found for the Qe (obtained from the above by 
setting the impact factor ha — 1), and then equation (|92|l is used to obtain the coefficient functions to any 
desired order in as- 

There is a slight caveat involving the resummed partons. In section 2, we defined the DIS(x) scheme at 
LL order, and gave the equivalent definition in the fixed order QCD perturbation expansion. However, this 
does not guarantee that the partons described in the two expansions arc the same. Note that this is not 
a problem at leading order in as, where all quantities are factorisation scheme independent. To be fully 
consistent with the conventional DIS scheme beyond LO, however, one must formulate the BFKL equation 
in 4 + 2e dimensions, regularising coUinear singularities via dimensional rcgularisation. It does not presently 
seem possible to achieve phenomenological results once the running coupling is introduced, although a recent 
discussion may be found in '47' . In the solution method used here, collinear factorisation is introduced in 
equation H84|) by modifying the lower limit of integration, with rcgularisation achieved by taking the gluon 
at the lower end of the BFKL ladder off-shell with virtuality Qq. This certainly does not correspond to 
the way in which the fixed order expansion is factorised. We believe, however, that the formal difference 
between the partons in the fixed order and small-x expansions can be neglected to a good approximation, 
and argue as follows. In the case of the LL BFKL equation with fixed coupling, solution by dimensional 
rcgularisation relates the gluon defined by solution of the BFKL equation (the so-called "Qo scheme" ) with 
the conventional MS scheme gluon |5II) : 

gQo ^ R^g^, (101) 

in A^-space, where Rn = 1 + 0{a'g). Thus the difference between the partons in the two expansions is irrel- 
evant when working up to 0{a'g) at fixed order. One hopes that running coupling effects do not seriously 
modify this situation. 

A recent discussion of scheme transformation between the Qq and MS schemes can be found in jlHli in which 
factorisation is derived in the dimensionally regularised BFKL equation (including the running coupling) 
using a saddle point approximation. The effect of the scheme transformation on a toy gluon distribution 
is investigated in j^HIi where it is noted that the size of the correction in transforming from the MS to Qq 
scheme is less significant than the difference between gluons obtained in the MS scheme by two rival groups 
(MRST and CTEQ). Bearing in mind that the transformation discussed in 0^1 is more involved than the 
shift from the Qo scheme to a DIS-like scheme considered here, one is justified in assuming the change in par- 
ton normalisation to be a small effect. One must also consider the quark sector, which is more complicated 
given that the effect of a transformation from the Qo scheme enters at the same order of the high energy 
expansion as the leading resummation in the splitting function Pqg. We assume, however, that as in the 
gluon case the effect of transforming from the Qp scheme DIS(x) quark to the quark distribution obtained 
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in the MS-subtracted DIS(x) scheme is small. 

There are other minor technical problems with the procedure specified here. Firstly, the FF resummed coef- 
ficient function is evaluated with (3o{nf), whereas the resummed VF coefficient and heavy matrix elements 
^Hg, are evaluated with /3o(?i/ + l). This results in a slight loss of continuity of i^2,ff across the matching 
point = M^. We checked, however, that this is an extremely small effect. Secondly, the coupling is im- 
plemented in the solution of the BFKL equation using the simple LO form as = l/(/3o log(fc^/A^)), whereas 
the implementation of the coupling when heavy flavours are accounted for is somewhat more complicated in 
that A changes discontinuously as one crosses each heavy quark threshold. This we also expect to be a small 
effect, and less significant than other effects due to choice of scale . 



3.2 The Resummed Heavy Flavour Coefficients 

To follow this procedure for the heavy flavour coefficients, one must expand the impact factors as a series 
in 7. However, they contain hypergeometric functions in (1 + AM^ / Q'^)~^ jJDj and so one must first Taylor 
expand them in a suitable Q^-dependent parameter. For low one can use z = {\ + AM"^ / Q'^)~^ and 
for high Q^, u ~ 1 — z. One needs sufficiently many orders of z or u as are needed to give an accurate 
description of the impact factor over the range of interest. The impact factors relevant to the FF coefficients 
are expanded in terms of z, and two expressions are needed for the VF factors for > A'P in terms of 
z and u respectively, overlapping at a suitable scale. We find that in the case of C|"j^, 18 orders of z are 
sufficient for a precision of a few parts permille up to = A/^. For C^g, 9 orders in z are needed up to 
z = 0.5 {Q^ = 4M^), and then 9 orders of u are sufficient for z > 0.5. For C^^'^^ , 8 orders of both z 
and u are necessary, again matching at z — 0.5. After these expansions, the resulting impact factors can be 
further expanded in 7. The series in 7 for J-'e must ultimately be truncated due to its being an asymptotic 
series. How many orders to include can be decided by examining the ratio ha/ha over the range 7 S (—1, 0), 
where ha is the double Taylor series in z (or u) and 7. Contributions from 7 < 1 are suppressed by inverse 
powers of the hard scale Q^. This ratio is required to be at least as accurate as the similar ratio of h factors 
for the massless case in I find that 5 orders of 7 are needed for each of the impact factors. Although 
the resulting structure functions are of finite order in as, the coefficients can be expanded to arbitrarily 
many orders in as/N. Following |24| . we require the series to be stable with respect to higher orders in as 
for N > 0.4 (corresponding to x' > 10~^). No significant improvement in any of the coefficients (to within 
0.1%) is found by taking more than 9 orders in as- 

The explicit expressions for the resummed A^-space functions are extremely lengthy (owing in part to the 
high number of terms in z or m needed when expanding the impact factors), and so are not reproduced here. 
They have the general form: 

9 n+1 3n—m, 

C^y-ias,N,u,) = E «SC„™(.)^, (102) 

n— m— 

with: 

5 

r s 

and w = z 01 u as appropriate. Note the presence of terms ^ ag/A^"+^. These are strictly absent in 
the full coefficients, but arise here from the truncation of the asymptotic series for the structure functions 
in double Mellin space and are weighted by inverse powers of /3o. This is justified in that the aim is 
to produce accurate phenomenological representations of the resummed coefficient functions with running 
coupling effects, and thus one must truncate the series as specified. The leading logarithmic terms are 

'^"See the discussion regarding the term fl^° of our equation 1901 in )24| . 
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Figure 4: Contour for numerical integration of equations and 



the same as those that would arise if running coupling corrections were not included, up to the order of 
truncation. The iV-space expressions are easily inverted to obtain the corresponding x'-space coefficient 
functions. For phenomenological applications, the coefficients C„m in equation (|102|l can be parameterised 
as a function of y = log(5^/4M^. We used: 

5 

Several ranges of y are needed to cover a range from the parton starting scale 1 GeV up to a TeV with 
an accuracy of less than a percent, but the resulting parameterised expressions are much quicker than the 
exact results when used in a global fit. 

No mention has yet been made of the power-suppressed correction to equation H9t)|l due to the ambiguities in 
its derivation. One can estimate this for each coefficient by comparing the analytically computed expression 
with the ratio of J-E.a and Qe &s found by numerical integration of equations H88|l and (|91|) in the 7-plane. 
Following , 24i . the contour for these numerical integrations is chosen as in figure 0| This is a deformation 
of the contour used in the definition of the inverse Mellin transform (parallel to the imaginary axis in the 7 
plane), and aids convergence of the numerical integration. One is allowed to make such a deformation given 
that branch points of the integrands of equations (|88l I91|l occur only along the negative real axis, which 
can be used to set the direction of all cuts. No singularities are then encountered in bending the inversion 
contour to that of figure 0] 

One can only perform the numerical integrations by choosing particular values of M and A. I use M = 1.5 
GeV and A = 150 MeV, where the heavy quark mass corresponds approximately to the charm quark. 
For each coefficient, the difference AC between the numerical and analytical results is found for N € 
{0.4,0.5,0.6,0.7,0.8,0.9,1,2,5,00} for two values t = to and t = ti, and then fitted to a function of the 
form: 

AC{N, t) = fit) exp i-At) + exp {~Bt) ^ 6„ (^77^) TTT-^- (^05) 



20 



0.3 



0.2 



0.1 




-5 -4 -3 -2 -I 

10 10 10 10 10 1 



Figure 5: The total resummed coefficient x' Cf^ [as,x' /M^) for t = 4.61 (Q^ = ~ 2.25GeV^) and 
Uf = 3 (solid line), together with the power-suppressed correction (dashed line) and perturbative piece 
(dot-dashed). Also shown is the LL correction with no running coupling corrections (dotted). 



The A'^-independent term is parameterised by finding AC at N — oo for various t values in the range 
to < t < ti. Usually a polynomial in t is sufficient for f{t). For C|^J we use to = 3.8 and ti — 4.61, 
corresponding at this A to = IGeV^ and = respectively. The A^-independent term was found 
using t £ {3.8, 3.85, 3.9, 4.0, 4.1, 4.2, 4.3, 4.4, 4.5, 4.61}, and the resulting expression for the x'-space correction 
is: 

x'ACiJ = ^ {(-1.5874- .94368 log(i- 3.6) -I- .33073t .587571og(i - 3.7))e~*^(l - x') 
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(106) 



with ^' = log(l/x'). Note that non-polynomial terms in the A^-independent contribution were needed for a 
reasonable fit with few parameters. A plot of the coefficient is shown in figureElfor t = 4.61, together with 
the power-suppressed correction and the result at LL order with no running coupling corrections. 



For CY^g, we used to = 4.61 and ti — 7, and the A-independent term was found by calculating the correction 
at t £ {4.61, 4.8, 5, 5.2, 5.4, 5.6, 5.8, 6, 6.2, 6.4, 6.6, 6.8, 7}. The power-suppressed correction in x'-space can be 
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Figure 6: The total resummed coefficient x' C'lg {as, x' , /IVP) for t = 4.61 {Q^ = ~ 2.25GeV^) and 
= 4 (solid line), together with the power-suppressed correction (dashed line) and perturbative piece 
(dot-dashed). Also shown is the LL result with no running coupling corrections (dotted). 



modelled by: 



x'A C^a = \ (-5.1401 + 4.3451t - mb2\t^ + 0.048160^^)6- 
■y 471" 



/.N \ 2.862916 

as(i) 



/ ,,\ \ 2.672627 / x 3.25314 / /.n x 3.811813 >/4 

262.40 ^'2^ 22.187 f^^^ C' - 7.1036 ^ "^^^^ ^ ^ 



Vas(4.61)y ■ Vas(4.61)y ' \as{4.M) J 4! 

/ X 4.02541 ^,5 / /,\ \ 4.24327 ^,6 

+3.9345 f^^^") -.58657 ^ "^^'^ ^ ^ 



Vas(4.61)y 5! Vas(4.61)y 6! 

and is shown with the variable flavour coefficient in figure 
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A plot of the VF coefficient as t increases is shown in figure Note that the power-suppressed correction 
remains a non-negligible correction to C^j^ even for intermediate x' values and higher t. This is due to the 
fact that both the correction and the coefficient are decaying as i — > oo, where both quantities are becoming 
less phenomenologically important. 

For the longitudinal coefficients, we used to = 4.5 and ti — Q (the correction dies away relative to the 
coefficient more rapidly in the longitudinal case, as Cj^^ does not vanish as Q'^/M'^ oo). The iV- 
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Figure 8: The x'-space coefficient x' C^^ for t = 4.5 and Uf = A (sofid fine), together with the power- 
suppressed correction (dashed) and perturbative piece (dot-dashed). Also shown is the LL result with no 
running couphng corrections (dotted). 



independent term was found by considering t e {4.5, 4.75, 5, 5.25, 5.5, 5.75, 6}. The result is: 



'ACL.g = ^ <i (-5.1401 + 4.3451t - .8652^^ 



47r 
1.3616 

.42836 



asjt) 
«5(4.5) 

asjt) 
as(4.61) 
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C - 1.9545 
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-3.72971 ,,4 
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(n.b. this is the correction to the FF coefficient for < A'P, and the VF coefficient for > M^). The 
coefficient and power-suppressed correction are shown in figure |S1 for t = 4.5. We note that the shape is 
very similar to C*^^ (figure E)) . There are also resummations contributing to the heavy quark coefficient 
Cl,hh — C'l where the form of the massless resummed quark coefficient can be found in pij . 



Once the gluon coefficients are known, the resummed coefficients to be convolved with the singlet distribution 
are calculated via the colour relation: 



(109) 



where Cp = 4/3, Ca = Nc = 3, and C'^\fg denotes the 0{as) contribution from the resummed gluon 
coefficient. This arises from the coupling of quarks to the BFKL 4-point function at the lower end, and is 



24 



true to LL order which is consistent with the order of the resummation. 

From figures one sees that the effect of the running coupiing is generaUy to suppress the high energy 
Hmit until lower values of x, and also to moderate the small x divergence. It is interesting to compare our 
results with those of alternative analytic approaches. In [SOI, the BFKL kernel is improved by requiring 
consistency of the small x expansion with the DGLAP limit of the gluon Green's function. A resummation 
of scale-dependent coUinear singularities of the kernel [SI] is implemented, and a generalised BFKL equation 
presented in which the kernels at each order of as depend on both N and 7 in double Mellin space. A 
method for solving this generalised equation (the "w-expansion" in the notation of (SUj) is given. This is 
further developed in [52, where results are presented for the gluon Green's function and gluon-gluon splitting 
function, as well as a discussion of impact factors needed for detailed phenomenology. A different approach 
is taken in |53| . where the BFKL anomalous dimension and kernel are improved using a duality relation 
between the two quantities derived on the assumption that the rightmost singularities of the gluon in the 
7- and N- planes coincide. The duality relation can be generalised to the running coupling case, and a 
proof that this can be extended to all orders was given recently |54j . Further improvements arise from the 
imposition of momentum conservation in the gluon evolution, and the resummation of coUinear poles |5T] 
- see |55| for a recent discussion. As far as phenomenology is concerned, however, there is no discussion of 
the incorporation of impact factors in the duality approach, which are needed for Pgg and the resummed 
coefficient functions. 

These approaches to the BFKL expansion are different to the method used here and in principle more 
complicated, but it is interesting to note that all three share the same qualitative results. The inclusion of 
the running coupling, present in all the approaches but most transparent in this one, stabilises the resummed 
perturbation theory and significantly moderates the small x divergence. Also common to each approach is 
the presence of a dip in the splitting function Pgg at intermediate x (see figure 6 of [S^ for the result adopting 
the method used here), before the small x limit sets in. Indeed, in all physical quantities the asymptotic 
small X behaviour seems to be pushed to lower values of x. 

3.3 The Heavy Flavour Matrix Element 

In order to implement the small x resummations in a fit to data, one also needs the heavy flavour matrix 
element Ang with running coupling corrections, evaluated at the matching scale. One thus follows the 
procedure above, using an "h- factor" given by equation H58(l with = M^. We again find the series in 
double Mellin space can be reliably truncated at 5 orders in 7, and the result for the resummed matrix 
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element is: 



" 47r 



.66667(5(1 ~x)- 1.0821(5s + (1.8609^ - 1.8609/3o) + (^-1.6174||- + 4.8521/3o^ 

-3.2347/3o) <5| + ^1.4987|j- - 11.594/3o||- + 24.290/3^^ - 14.195/3;^^ <5| + (^7.8212|j 

-51.369/3o^ + 121.63/3^^ - 122.63/3;^^ + 44.549/3;^ + -1.7449^ + 16.204^ 
3! 2! / \ /5o6! 5! 

-61.577;3o^ + 118.82^2^ - 114.39^^^ + 42.685^0^^^ a% + (^1.9728i^ - 26.10l|^ 
+135.73/3o|t - 357.64/3^^ + 503.70/3i^|^ - 358.74/3^^ + 101.07/?^^ "5 + 5.5968- 



.61.77l|^ + 282.34/3o|J - 687.81/32|^ + 955.67^^^! _ 748.174^1! ^ 300.36^^^ 



-46.201;3o'O 4 + (^-1-3983^ + 20.599^ - 138.35|j + 564.42^o|y - 1560.3/32|^ 
+3036.6/33 1! _ 4111.0^4^ ^ 3652.2^5^ _ i886.7^o'f + 423.92^0^^") 
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where ^ = log(l/x). The power-suppressed correction is somewhat easier for this quantity, as it is only 
needed for one value of t. The difference between the analytical and numerically evaluated results is 
found for N £ {0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1, 2, 5, oo} and fitted to a function of N only. The result in x-space 
is then: 

£2 £3 M ^5 

xA^ = -0.013550(5(1 -x)- 0.058138 + .39486^ - .62386^ + .41538^ - .12522-^ + 0.014138^. (Ill) 

2! 3! 4! 5! 

A plot of the matrix element is shown in figure |3 alongside the corresponding LL order result evaluated 
to 0{ai.). Note that the power-suppressed correction dominates the analytic result. This problem is made 
worse if one includes more powers of 7 in double Mellin space due to the instability of the asymptotic 
expansion. Neither is the size of the correction reduced by expanding the analytic result to higher powers 
in as- However, this is not really a problem given that the power-suppressed correction is modelled very 
well due to only being required at a single value of Q^. One can see from figure El that the effect of the 
correction is to restore positivity of the matrix element at small x, as noted in iV-space previously. We again 
note that the running coupling corrections moderate the small a; divergence from the strictly LL result. The 
resummed singlet matrix element Anq is calculated from Ahq analogously to equation p09|) . 



4 Global fit to Scattering Data 
4.1 Data Sets and Details of the Fit 

In this section we present the results of a global parton fit adopting the DIS(x) factorisation scheme outlined 
in the preceding sections, implementing the resummed splitting functions in the parton evolution and the 
resummed coefficient functions in the theoretical F2 and predictions for electromagnetic current DIS. We 
include structure function data from the HI [T1[T31 US] and ZEUS [HI HOI ES| collaborations at HERA; 
proton data from BCDMS NMC SLAG and E665 ED; deuterium data from BCDMS iSl, 

NMC, SLAG and E665; GGFR data on F:^';^^ {x,Q^) gSlEl; data on the deuterium-proton ratio /F| 
from NMG [HSl; charged current data from HI Jl] and ZEUS data on the charm structure function 
F2.C from HI m and ZEUS [TT]. The non-DIS data sets used are Drell-Yan (DY) data from the E866/NuSea 
collaboration DY asymmetry data from NA51 data on the DY ratio 

(T^y/agV from E866 EH; 
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Figure 9: The total resummed matrix element Ang with running coupling corrections for n/ = 4 (solid line), 
together with the power-suppressed correction (dashed) and perturbative piece (dot-dashed). Also shown is 
the LL result with no /Jo-dependent corrections (dotted). 
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W-asymmetry data from CDF [TOj. By means of a consistency check, we also compare the resulting gluons 
with a gluon from a previous NLO fit, which includes Tevatron jet data from the DO "717 and CDF 
collaborations. 

The parton distributions are parameterised using the forms given in |73| at a starting scale of Q|„ji = IGeV^, 
allowing for a valence-like or even negative gluon at the starting scale. The running coupling is implemented 
according to the prescription of 25 for dealing with heavy flavour thresholds. As well as including small x 
resummations as outlined in this paper, we also include a resummation of the leading large x divergences 
in the anomalous dimension Pqq by modifying the argument of as acting on the plus distribution from 
—^ Q^{1 — x) |75j. As will be seen, such a resummation is needed in order to achieve agreement between 
theory and data at high x. To avoid encountering the Landau pole in the coupling as a result of this modi- 
fication, as is frozen for < 0.144GeV^ (as ~ 0.6). 

A problem occurs with the Drell-Yan data in that the NLO QCD corrections to the cross-section are known 
to be large [2^1 • To obtain a reasonable fit to the Drell-Yan data, we thus employ a constant "K- factor" 
(~ 1.4) as in [77], whose effect is normalise the DY data outside the confines of the systematic error. Al- 
though strictly absent at LO, such a factor does not have any xi, x^ or dependence and thus only poorly 
models the known higher order corrections. 

4.2 Results 

In table 121 we show the values for each data set obtained in the LL fit, together with the results obtained 
from a LO fit with no resummations at small or large x. 

The LO fit is reasonable, but misleading. It fails due to the fact that the evolution is too slow. Some 
compensation is achieved by increasing the QCD scale parameter A, in that by raising the coupling constant 
as the quarks evolve faster, and one can increase the theoretical predictions to be more in line with the 
data. The LO fit gives as(Mf ) = 0.1305, which is indeed rather high when compared to the world average 
of 0.1187(20) [7S|. This is not strictly a problem at LO, where one is free to redefine the renormalisation 
scale - chosen here to be = ■ It is still, however, a cause for concern given that is the natural scale 
choice. As well as a raised coupling constant, the HI and ZEUS normalisations are required to be lower in 
the LO fit than the corresponding resummed fit parameters. The description of F2 at small and high x fails 
ultimately because the shape is too flat. The fit to the SLAC data, for example, improves significantly once 
high X resummations are included. 

The description of the small x DIS data is improved by resummation. In figures ^1 and ^] we show the 
resummed theoretical predictions alongside the data for a range of low x values. The data are clearly fit very 
well, with no systematic tendency to undershoot. This is in contrast to results from a NLO fit, also shown 
in figures [Till and ITTI The LL fit clearly has the increased slope needed to continue to fit the data well as 
increases. The charm structure function F2 is also fit well, and in figure [T^ we show the resummed theoret- 
ical prediction alongside the data and NLO fit results. Both the LL and NLO fits give a good overall fit, due 
to the somewhat large uncertainties associated with many of the data points. However, the resummed fit 
performs better at very small x. Even so, theory still underestimates F| at the lowest x values (< 3 x 10^"*). 
This may indicate the importance of higher orders in the resummed expansion. 

Note that a poor fit is obtained for the DY data in the LL fit. This is not surprising, given the large 
perturbative corrections at NLO and the crude nature of the constant factor K.^^ . Consistent with this 
expectation, the NLO fit performs much better for this data set. The DY ratio data is better fit, which one 
expects given that some perturbative uncertainty cancels in taking the ratio of cross-sections. 
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Figure 10: Theoretical predictions for the structure function F2 alongside the data (normalisation dictated 
by the fit), for 5 x IQ-^ < a; < 4 x 10"*. 
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Figure 11: Theoretical predictions for the structure function F2 alongside the data (normalisation dictated 
by the fit), for 5 x 10"'' < a; < 8 x IQ-^. 
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Figure 12: Resummed predictions for the charm structure function F2 alongside HERA data, for 1.3 x 10 < 
x<2x 10"^ 
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Data Set 


No. data pts 




xio 


Xnlo 


HI ep 


417 


342 


414 


427 


ZEUS ep 


356 


282 


287 


279 




27 


26 


24 


32 


BCDMS n p 


167 


170 


263 


191 


BCDMS fi D 


155 


230 


208 


216 


NMC ^ p 


126 


111 


154 


136 


NMC n D 


126 


89 


141 


103 


SLAC ep 


53 


77 


195 


67 


SLAC eD 


54 


74 


188 


56 


E665 /i p 


59 


59 


56 


61 


E665 M D 


57 


55 


54 


51 


CCFR F^'^ 


74 


99 


158 


83 


CCFR if ^ 


105 


138 


126 


115 


HI CC 


28 


33 


34 


29 


ZEUS CC 


30 


46 


34 


35 


NMC n/p 


156 


164 


144 


154 


E866/ NuSea DY 


174 


307 


276 


237 


NA51 DY asym. 


1 


11 


6 


11 


E866 afy/aPPy 


15 


7 


29 


10 


CDF W asym. 


11 


16 


26 


14 


Total 


2181 


2336 


2817 


2307 



Table 2: The quality of fit from the LL and LO fits for each dataset, as well as results from a previous NLO 



Whilst the small x fits to the HERA data improve, one can see from table |21 that fits to some of the data 
sets actually worsen when resummations are implemented. One expects some tension between theory and 
data for non-electromagnetic DIS and other processes, due to the fact that resummed partons are being used 
without the appropriate resummed impact factors. However, these data sets do not contain points at very 
low values of x. Instead the problem lies at intermediate values of x (10~^ ^ a; < 5 x 10~^), where the 
theoretical prediction underestimates the deuterium data (proton structure function data, on the other hand, 
is fit rather well - in many cases the LL fit to F2 outperforms the NLO fit). This indicates that the effect of 
the resummation in the partons is felt at higher values of x than one naively expects. The higher evolution in 
the moderate x region means that the sea quarks evolve much more quickly than the non-singlet quark com- 
binations, leading to problems describing the relative shape of the sea and valence quarks. This accounts for 
the poorer fits to the deuterium and charged current data, and also the results for the Drell-Yan data. One 
sees from table|21that the E866 DY set is fit worse in the LL fit, even with the use of a variable K-factor. The 
exaggerated influence of small x resummations is also evident in the value of as{M'^) obtained in the LL flt of 
0.1126, which is rather low compared to the world average. For comparison, the NLO fit gives as — 0.120 

We find that the gluon obtained from the LL fit is inconsistent with the gluon obtained from a previous 
NLO fit to the Tevatron jet data. The LL gluon is shown in figure IT^ together with a typical NLO gluon 
for two different scales. At the starting scale in a LO fit, the gluon remains positive at small x, due to the 
fact that larger partons together with a raised coupling constant are needed to try and fit the HERA data. 
Once higher order corrections or resummations are implemented, one expects a lower gluon at small x, due 
to the faster evolution in this regime. In fact, figure IT!^ shows that a negative gluon is needed at the starting 
scale in both the LL and NLO fits. However, the LL gluon turns negative at significantly higher x (different 
by about an order of magnitude) than the NLO gluon. Furthermore, there is a large suppression of the LL 

^^One must boar in mind, however, that the definition of Qg is order dependent. 
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Figure 13: The gluon distribution at Qf^u = IGeV^ and — lOOGeV^ obtained in the LL resummed fit 
and a previous NLO fit [JHI- Also shown is the result obtained from the modified LL fit described in the 
text. 
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result relative to the NLO gluon in the intermediate x region. This is further evidence that the effects of 
the small x resummations are felt too strongly at moderate x. The suppression of the gluon in the LL fit 
persists at higher Q^, as can be seen from the right-hand panel of figure [TSl We note that the left-hand plot 
shows the input scale gluon from the resummed fit agreeing extremely closely with the corresponding NLO 
gluon at small x. This is nothing more than a coincidence, as can be seen from the right-hand panel where 
the evolution drives the gluons to be quite different at higher scales. 

The solution to the problem of resummations manifest at moderate x is found by considering higher order 
terms in the fixed order expansion - NLO corrections help to suppress small x resummation effects until 
lower values of x, due to the inclusion of the correct higher order moderate x behaviour. Indeed, a NLO 
fit with no resummations shows significant improvement over a LO fit, and the subsequent improvement 
upon performing a NNLO analysis stems largely from the inclusion of the extra small and large x terms 
(see [Zni)- The characteristic feature of the complete NLO Pqg and Pqq is a dip at moderate x, followed by 
a rapid increase at very high x due to the presence of high x divergences (absent in the MS scheme). One 
can therefore approximate subleading behaviour at small x by modifying the resummed splitting function, 
replacing the 0{ag) term as follows: 



^LL,mod.{2) ^ 34 



f A 

x a;" 



(112) 



where one chooses the values of A and a to closely model the dip at moderate x. We use a — 0.47 and 
A = 1.2. Note that this does not have a serious impact on the splitting function at high x, given: 

lim(a;"i - 1.2a;-") = -0.2, (113) 

a:— >1 

whereas the complete NLO DIS scheme splitting function is divergent as a; ^ 1. 



We performed a global fit using this modified anomalous dimension (which also affects Pqq by equation 
(|109|l ). No modification to the resummed Pgg (and, thus, Pgq by the colour relation) was implemented. 
However, this has no term at 0(q;|) in the resummed expression. Furthermore, the NLO corrections in the 
fixed order quantity are small. We find a more favourable comparison to the gluon arising from the Tevatron 
jet data. The gluon which emerges is shown alongside those from the other fits in figure 1131 and one can 
see that although it is qualitatively the same as the LL result, the suppression at moderate x is no longer 
as severe. Indeed, it is qualitatively more similar to the gluon obtained from a NLO or NNLO fit. This is 
therefore an indication that a next-to-leading order analysis (n.b. NLO and NLL) is what is really needed. 
It is still not possible even in the modified fit to obtain a very good fit to all the charged current data (or, 
for that matter, the gluon from the jet data), suggesting again the necessary inclusion of higher orders. 



5 The Longitudinal Structure Function 

In this section we discuss the resummed predictions for ■ A plot of the value of the longitudinal structure 
function obtained from the LL resummed fit described in the previous section is shown for various values of 
in figure El Also shown are the results obtained from the LO fit undertaken in this paper, together with 
similar results from previous NLO and NNLO fits [50] . One sees that the LO estimate is far higher at small 
X when compared with higher orders in the fixed order perturbation theory, owing to the larger gluon at 
low and the higher value of as needed to fit the F2 data. Comparison with the NLO and NNLO curves 
reveals the apparent perturbative instability at small x. This is particularly noticeable in the = 2GeV 
plot, corresponding to a scale at which the gluon is negative (at NLO and NNLO) at low x. One needs the 
divergent low x terms at 0{a^) in the longitudinal coefficient functions to counteract the negative small x 
gluon leading to positivity of the structure function (see SHI for a discussion). The poor convergence of the 
fixed order calculations at low x (of greater significance at low Q^) indicates that high energy resummations 
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Figure 15: Resummed prediction for the longitudinal structure function as varies, for a: = 10 ^ and 
X = 10"''. Also shown are fixed order results, and the result from the dipole fit of |5T]. 



are probably necessary in this regime. 

One sees that at small x and Q^, the resummed result lies above the NNLO fit, but evolves more slowly than 
the fixed order fits as increases. The important point to note is that it is much more stable than the fixed 
order results at small x, and lies somewhere between the LO and NLO results. A less desirable feature of 
the resummed fit, however, is the behaviour at moderate and high x. In these regimes the resummed result 
is severely below the fixed order results, even as reaches values at which the NLO and NNLO results are 
closer together (n.b. denoting less sensitivity to the small x instability). The reason for this, as noted in 
the previous section, is that the gluon from the resummed fit is much smaller at higher x than should really 
be the case. Also, the lower value of A generated by the resummed fit decreases the contributions from the 
quarks at higher x. 

To emphasise the improvement at low x values in the resummed fit, we show in figure ^1 the evolution of Fi^ 
for X = 10~^ and x=10^'*. The instability in the fixed order results is clearly visible at lower Q^, where the 
resummed fit is much more sensible. The resummed result is much fiatter in than the fixed order results. 

It is interesting to compare the resummed fit performed here where the evolution is set up in terms of parton 
densities, to a previous resummed fit which evolves the structure functions directly |17| using physical 
anomalous dimensions ¥F and which had a significantly less rigorous treatment of heavy flavours. From 
figure El one sees that the results are qualitatively consistent between the two fits. They converge at very 
high a; to a value lower than the fixed order results, given that they both suffer from a gluon which is too 
small in this region. Also, NLO and NNLO corrections to the longitudinal coefficient fimctions are important 
at high X, and these are missing from the resummed fit. However, the evolution is moderated a little at 
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intermediate x in the physical anomalous dimension fit, due to an additional convolution in evaluating the 
derivative of the structure functions. For example, one has for the singlet structure function: 



dF2{x,Q^) 
91ogQ2 



T2l{x) ® Fl{x, Q^) + T22{x) F2{X, 

^2l{x)(E) [CLg<E)g{x,Q^) + CLq<»^{x,Q^)] + ... 



(114) 



(115) 



where the ellipsis denotes the term in i^2- The term in the gluon distribution thus involves an intermediate 
convolution with the coefficient: 



This tends to zero as x — *■ 1, and the effect of this additional convolution is to suppress the phenomenological 
effect of the small x resummations in the physical anomalous dimension T2L to lower values of x. The result 
is that the effective gluon and in such a fit are increased in the moderate x region of x ~ 10^^ — 10^^. 
Both the resummed fits are lower than the fixed order results at moderate x, as can be expected from the 
suppressed gluon. 

In figure ^1 we also show the value of Fl obtained from the dipole model fit presented in , itself 
similar in approach to previous dipole fits |82[ 1831 IMl I85| . There, the structure functions are written as 
the convolution of a dipole cross-section (dependent on the unintegrated gluon distribution) and a dipole 
wavefunction representing the probability of a photon fluctuating into a quark-antiquark pair. This is 
essentially a restatement of the LL kr factorisation theorem of equation (|23|l . although the gluon distribution 
is modelled differently using a parameterisation based on evolution according to the approximate anomalous 
dimension: 



One sees from the figure that at low the dipole fit is very close to the LL resummed fit performed 
here. They both share the feature of a gluon which is much too small, compared to a DGLAP gluon, at 
moderate and high x - where the DGLAP results should be reliable. As increases, however, the dipole 
result overshoots the resummed fit, owing to its much higher evolution. This can also be seen in figure 
1151 where the behaviour of the dipole result as increases is shown. It must be stressed that both the 
dipole and LL resummed fits give very good descriptions of the low x data. The shortcoming of both, 
however, is that they are not quantitatively correct over the whole x range, owing to the need to suppress 
the small x resummation effects at higher x. Note that the kx factorisation approach is easier to extend to 
higher orders. At LL order it is formally equivalent to the dipole approach. Beyond LL order, the dipole 
approach suffers conceptual problems. For example, including the exact gluon kinematics in the photon- 
gluon impact factor (a significant effect at NLL and higher orders) leads to the non-conservation of impact 
parameters describing transverse dipole size throughout the interaction |86j . which is a necessary property 
for the applicability of dipole arguments. Although it is possible this property may be recovered by inclusion 
of full NLL effects in the dipole wavefunction, there are no such conceptual barriers to extending the kx 
factorisation approach. Furthermore, there are no difficulties in matching to the fixed order expansion in 
the factorisation framework. It is also more suited to the consideration of proton-proton collisions, which do 
not have a direct interpretation in terms of dipoles. 

6 Conclusions 

In this paper we have proposed a new variable fiavour number scheme, the DIS(x) scheme, for use in deep 
inelastic scattering. The motivation for this scheme is to provide a way of carrying out heavy flavour calcu- 
lations at small a; in a form that can be consistently implemented alongside the flxed order QCD expansion. 
We have shown how to obtain the resummed coefficients and heavy matrix elements at small x in double 



CLg{x)^-^ix{l-x)+0{al). 



(116) 




(117) 
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Figure 16: Resummed prediction for alongside the result obtained from a previous resummed fit using 
physical anomalous dimensions. Also shown is the structure function obtained from the dipole model fit of 
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Mellin space from existing expressions for the photon impact factors at LL order. One then obtains results 
in X and space after solving the BFKL equation. Here we have adopted a particular solution method 
which resums running coupling corrections, but the formulation of the DIS(x) scheme in {N, 7) space is more 
general than this and can be used with other BFKL approaches at small x. Alternatively, it can be used at 
fixed order even if one is not considering high energy resummation. 

The advantage of the method adopted here is that it produces analytic results which can be directly used 
in parton fitting routines. Furthermore, /3o-dependent corrections are effectively resummed to all orders (at 
least, to any desired phenomenological precision). It offers a more physically motivated definition of both 
the massless partons (due to the use of the DIS scheme) , and the heavy partons (due to the choice of heavy 
quark coefficients). We used our scheme in a LO global fit to scattering data, including LL resummations 
with running coupling corrections. The running coupling is seen to significantly moderate the strong high 
energy divergence observed in physical quantities in a purely LL fit. A good overall fit is obtained, and the 
fit to HERA data at small x is improved relative to the LO fit with no resummations. The description of 
the HERA charm data for F2 is also improved. Compared with a NLO fit, the LL results give a significantly 
better description of the low x data and describe most of the proton data well. However, the fit worsens 
at moderate and high x for the non-proton data, due to the fact that the small x resummations lead to a 
decrease of the coupling constant, and an incorrect relative shape between the sea and valence quarks due to 
the increased evolution of the former. There is also a very small gluon distribution in the moderate x region, 
which is inconsistent with previous gluons obtained from Tevatron jet data. This suggests that higher order 
corrections in both the fixed order and resummed expansions are what is really needed to understand the 
data. The effect of higher order corrections is to suppress the influence of small x resummations to lower 
values of x. We demonstrated this here by making a small modification to the LL fit, which tempers the 
increased evolution of the quarks at moderate x. This resulted in a gluon distribution more in line with that 
obtained from the Tevatron jet data. The gluon from the modified fit shows a significant increase in the 
moderate x region, as expected, which allows for a better description of the high x data. 

One does of course expect some degree of tension in the resummed fit between the theoretical predictions 
and the data for non-DIS sets. This is due to the use of partons obtained using small x resummed evolution 
kernels, but without including resummed coefficients for e.g. charged current scattering. This should not 
be too significant in practice, however, given that the non-DIS data sets do not include data points at very 
low X. The fact that the non-DIS data are not well fit is itself an indication that NLO corrections are needed. 

Our resummed results for the longitudinal structure function were presented in section 5. Importantly, the 
resummed fit produces much more stable results at low x than the fixed order results, particularly at low 
where resummations are expected to be more important due to the higher coupling constant. However, 
at moderate and high x the resummed results undershoot the fixed order results. This feature is shared at 
lower by the dipole model fit of which together with the LL fit leads to a gluon which is much too 
small at moderate x. A previous resummed fit using physical anomalous dimensions gives similar qualitative 
results. One does not expect to match the NLO and NNLO results at high x in the LL fit, due to missing 
high X terms in the longitudinal coefficients. The results at low x are very encouraging, however, and give 
a strong indication that high energy resummations are a necessary prerequisite for a sensible prediction of Fj^ . 

Our final conclusion is therefore that, although we have demonstrated the need for high energy resummations 
in order to correctly describe structure function data at small x, higher order corrections in the fixed order 
expansion are needed to achieve a truly quantitative agreement over the full x range. For consistency one 
should build NLL resummations over a NLO fixed order expansion. Our investigation into implementing 
such a scheme in an approximate framework is ongoing [S7| . 
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